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We investigated the magnetism under a magnetic field in tlie quasi-two-dimensional organic Mott 
insulator K-(BEDT-TTF)2Cu[N(CN)2lCl through magnetization and "C-NMR measurements. We 
found that in the nominally paramagnetic phase (i.e., above Neel temperature) the field-induced 
local moments have a staggered component perpendicular to the applied field. As a result, the anti- 
ferromagnetic transition well defined at a zero field becomes crossover under a finite field. This un- 
conventional behavior is qualitatively reproduced by the molecular-field calculation for Hamiltonian 
including the exchange, Dzyaloshinsky-Moriya (DM), and Zeeman interactions. This calculation 
also explains other unconventional magnetic features in «;-(BEDT-TTF)2Cu[N(CN)2]Cl reported in 
the literature. The present results highlight the importance of the DM interaction in field-induced 
magnetism in a nominally paramagnetic phase, especially in low-dimensional spin systems. 



PACS numbers; 75.30.Kz, 74.70.Kn, 74.25.Nf, 68.35.Rh 



I. INTRODUCTION 

The quasi-two-dimensional (quasi-2D) organic con- 
ductor, K-(BEDT-TTF)2X, exhibits various phases, de- 
pending on anion X, temperature and pressure, and 
have provided a good experimental stage for inves- 
tigating fundamental problems in condensed-matter 
physics^ [BEDT-TTF, abbreviated as ET hereafter, is 
bis(ethylenedithio)tctrathiafulvalene]. For instance, n- 
(ET)2Cu[N(CN)2]Br and k-(ET)2Cu(NCS)2 show un- 
conventional superconductivity at ambient pressurCf?'- 
and intensive studies regarding its origin are in progress.^ 
k-(ET)2Cu[N(CN)2]C1 (hereafter abbreviated as k-C1) 
is a Mott insulator and suitable for the study of the 
bandwidth-controlled Mott transition, because this ma- 
terial undergoes the Mott transition by soft pressure 
(~ 25 MPa)i^ In fact, recent experiments on k-C1 us- 
ing the pressure-sweep technique have revealed funda- 
mental aspects of the Mott transition: the first-order 
transition with a finite-temperature critical endpoint 
(Tcr ~ 40 K)^iii2ii2ii£ and unconventional critical expo- 
nents of the Mott criticalityi^ Another interesting as- 
pect in K-(ET)2X is the magnetism on the frustrated 
triangular lattice. As will be seen in Sec. II, n- 
(ET)2X has an anisotropic triangular lattice, where lo- 
calized moments with antiferromagnetic (AF) correla- 
tion are subject to the so-called spin frustration. In 
this context, the Mott insulator k-(ET)2Cu2(CN)3 is 
an intriguing material, because it has nearly isotropic 
triangular lattice. In fact, neither the Neel order nor 
spin-gapped behavior are observed down to 32 mK in 
NMR^^'^'^ and /iSR measurements A quantum spin 
liquid has been suggested as the possible ground state 
of K-(ET)2Cu2(CN)3» 

In this paper, we focus on the magnetism of the Mott 
insulator, k-C1. Now, its magnetism is understood as 
follows: the AF transition occurs at Tn 27 K accom- 
panied by weak ferromagnetism, which is due to the cant- 
ing of antiferromagnetically ordered spins;^ This canting 



is attributed to the Dzyaloshinsky-Moriya (DM) interac- 
tion, which is inherent in the k-(ET)2X (see Sec. II). 
However, further investigations seem to be needed. For 
instance, Hamad et al. found two curious features in 
^H-NMR measurements:^® with increasing magnetic field 
from 0.58 T to 9.3 T, the peak structure of spin-lattice 
relaxation rate 1/Ti around Tn is strongly suppressed, 
and the peak temperature (usually regarded as Tn) in- 
creases from ^ 22 K to ~ 28 K, although in conventional 
antiferromagnets 1/Ti is field-insensitive and Tn is ro- 
bust (or decreases under strong field). Hamad et ali^ 
interpreted the field-dependent 1/Ti as a consequence of 
slow spin dynamics due to the spin frustration; however, 
the unexpected increase in the nominal Tn with magnetic 
field remains puzzling. 



To get more insight into the magnetism of k-CI, we 
combined magnetization and NMR measurements. On 
the basis of the experimental results, we argue that 
the AF transition at a zero magnetic field changes into 
crossover under a finite field. This statement leads to an 
intuitively strange consequence that field-induced local 
spin moments in the paramagnetic phase have a stag- 
gered component as well as the uniform one, but we 
actually found the non-coUinear state under a magnetic 
field. We also performed the molecular-field analysis of 
Heisenberg model including the DM and Zeeman inter- 
actions and succeeded in reproducing qualitatively not 
only our findings but also the two curious features found 
by Hamad et alJ^ This means that the above-mentioned 
puzzling magnetism under magnetic field is the manifes- 
tation of the interplay between DM and Zeeman inter- 
actions. To our knowledge, the importance of DM inter- 
action in a nominally paramagnetic phase has not been 
recognized so far. In this paper, we demonstrate that 
the DM interaction plays a key role in the magnetism 
under magnetic field even above Tn, especially in low- 
dimensional spin systems. 
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FIG. 1: (Color online) (a) Side view of crystal structure of 
K-Cl. (b) Top view of schematic crystal structure of the con- 
ducting ET layer. The pair of lines represents a (ET)2 dimer. 
(c) Dimer-lattice model of the conducting layer in k-C1. (d) 
Localized spin model applicable to k-C1. 



II. CRYSTAL STRUCTURE OF k-(ET)2X 

«;-(ET)2X has a quasi-2D layered structure composed 
of conducting ET layers and insulating anion layers [Fig. 
[IJa)]. In the conducting layer, ET molecules form dimers 
with face-to-face configuration [Fig. [Ijb)]. The mono- 
valent anion X introduces a hole into an antibonding 
dimer orbital, which forms a two-dimensional half-filled 
band through the inter-dimer transfer integrals. Thus 
the dimer lattice can be reduced to an isosceles triangu- 
lar lattice as shown in Fig. [TJc), which is characterized by 
two inter-dimer transfer integrals, t and t' . In the case 
of the Mott insulating k-(ET)2X, such as k-C1 (space- 
group Pnma), a hole is localized at a dimer. Because of 
the triangular lattice, the localized spins with AF correla- 
tion are subject to the spin frustration, of which strength 
is characterized by the ratio t' /t. The t' /t of k-C1 is es- 
timated at ~ 0.75 (i.e., anisotropic triangular lattice)ii^ 
The unit cell of k-C1 extends over two adjacent conduct- 
ing layers [layer A and layer B, see Fig. [TJa)] and the 
each layer includes two inequivalent dimers [dimer 1 and 
dimer 2, see Fig. [IJb)]. Thus, the unit cell turns out to 
contain four inequivalent dimers: Al and A2 in layer A, 
and Bl and B2 in layer B. These notations are after Ref. 
19. 

When the magnetism of k-C1 is discussed in terms of 
the localized spin model, t and t' are replaced by AF 
exchange interactions, J and J', respectively. However, 
one should keep in mind that the DM interaction also 
works between the nearest dimers [Fig. (Hd)] because 
the local inversion symmetry is absent between dimer 
1 and dimer 2 [see Fig. [Ub)]. In fact, as seen below, 
the puzzling magnetism found under a magnetic field is 
understood as a consequence of the DM interaction. 



III. EXPERIMENTAL 

A. Magnetization measurements 

To know the macroscopic magnetism of k-C1, we per- 
formed magnetization measurements at and 7 T with 
a SQUID magnetometer. In the measurement of spon- 
taneous weak ferromagnetism under T, k-C1 was first 
cooled down to 2 K under 0.1 T (parallel to the con- 
ducting layer) and then the magnetic field was switched 
off. Finally, the magnetization along the plane was mea- 
sured in an ascending temperature process under T. 
The uniform spin susceptibility Xspin at 7 T was obtained 
after subtracting the core diamagnetism [—4.7x10^'* 
emu/(mol f.u.)] from the measured susceptibility, which 
was determined by the magnetization divided by the ap- 
plied magnetic field. 

B. NMR measurements 

In the NMR measurements, we used a k-C\ crystal 
where two central carbon sites in ET molecule are sub- 
stituted by ^•^C isotope (hereafter abbreviated as k-^^C2- 
Cl). Under the field applied in an arbitrary direction, k- 
^^C2-C1 shows 16 ^^C resonance lines at most. This line 
profile arises from the following three factors.^*' First, as 
described above, there are four inequivalent dimers in the 
unit cell: Al, A2, Bl, and B2. Second, the inclined face- 
to-face geometry of ET molecules in the dimer makes the 
two central -'^^C sites in ET inequivalent [Fig. [H^a)]: the 
^^C site closer to the center of the dimer is labeled with 
"inner", and the other ^'^C is labeled with "outer". Fi- 
nally, a nuclear dipolar splitting between the central ^^C 
sites doubles the number of lines. Below we deal mainly 
with the outer site, since the NMR properties of the inner 
and outer sites show the similar behavior. 

In the present NMR measurements, the magnetic field 
H of 7.4 T was applied parallel to the a axis. In that 
case, the NMR shift of the outer site at dimer i (= Al, 
A2, Bl, B2), (in ppm), is given by 

j\fa J\,fb i\/fc 

Kf = + + Ar-^ + , (1) 

where M" denotes the a-axis (a — a, 6, c) component of 
the local spin moment Mi at dimer i, A^" is a compo- 
nent of the hyperfine coupling tensor of the outer site in 
dimer i, and Cf (in ppm) is the chemical shift of the site. 
As discussed by Smith et al,^S, A'^ is positive for every 
a, and the hyperfine coupling constants at other dimers 
can be determined by applying the symmetry operation 
inherent in k-C1 to A'^ (space-group Pnma). The de- 
rived relations between the hyperfine coupling constants 
at different dimers are as follows: 

^Ai = <2 = -A%\ - -<2 > 0, (3) 
AZ = -A'X2 = A%\ = -A%^2 > 0. (4) 
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FIG. 2: (Color online) (a) Side view of ET dimer. The cross 
mark represents the center of the dimer. (b) ^"^C-NMR spec- 
trum of K-Cl under H\\a at room temperature. 



Under J?||a, the dipolar splitting nearly vanishes be- 
cause H approximately forms the so-called magic angle 
(~ 54.7°) against the 13C="C vector in ET by hap- 
penstance. At room temperature (in the paramagnetic 
state), all dimers are equivalent under H\\a] thus ^^C- 
NMR spectra of k-^^C2-C1 shows only two lines coming 
from the inner ^^C and outer ^'^C [Fig. [DJb)]. A marked 
feature in the experimental results is that the two-line 
shape holds even at low temperatures where the spins are 
ordered. This means that Kj^-^ — K'^^ — — ^B2 
satisfied in the whole temperature range. By using Eqs. 
(l)-(4), this relation is rewritten as follows: 



Mil = = Af^i - 
Ml, = -Ml, = M%, = -All,,. 



(5) 
(6) 
(7) 



Equations (5)-(7) put strict constraints on the spin con- 
figuration under H\\a. In Sec. IV, we determine the spin 
structure at low temperatures on the basis of these con- 
straints. The absolute values of the Knight shift in the 
present study are calibrated with the room-temperature 
line position of the outer site, which is determined to be 
336.2 ppm by Smith et alJ^ 



IV. EXPERIMENTAL RESULTS 

In this section, by investigating the macroscopic and 
microscopic magnetisms of /c-Cl, we demonstrate (i) the 
increasing "Tn" with a magnetic-field and (ii) the field- 
induced non-coUinear local moment in a paramagnetic 
phase. First we show the macroscopic magnetism. The 
temperature dependence of Xspin of a k-C1 single crys- 
tal measured at 7 T (parallel to the conducting layer) is 
presented in the inset of Fig. [31 The result is consistent 
with Xspin of the powdered sample measured at 0.1 T.^^ 
The weak ferromagnetism observed below 30 K is due to 
the DM interaction and therefore should appear sponta- 
neously even under a zero field. In fact, as shown in the 
main panel of Fig. [3l the magnetization measurements at 
a zero field confirm a spontaneous weak ferromagnetism 
parallel to the conducting plane below ~ 23 K. This tem- 
perature is regarded as Tn at a zero field (abbreviated 
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FIG. 3: Temperature dependence of magnetization under a 
zero magnetic field. Inset: Temperature dependence of spin 
susceptibility measured under 7 T (parallel to the conducting 
layer) . 



as the zero-field Tn hereafter) and is in good agreement 
with the peak temperature of 1/Ti under 0.58 T.^- Note 
that the zero-field Tn is appreciably lower than the peak 
temperature under 3.7 T (~ 27 K)M- Thus, the present 
results support the finding by Hamad et alM that "Tn" 
seems to increase with a magnetic field (actually, as seen 
below, the peak temperature of 1/Ti should be regarded 
as crossover temperature rather than ordering tempera- 
ture in K-Cl). 

Next, we turn to the microscopic magnetism of k-CI. 
Figure HJa) shows the temperature dependence of ^^C- 
NMR shift at the outer site. The large shift toward low 
temperatures refiects the growth of local moment, con- 
sistent with the previous report li^^, As demonstrated by 
Smith et al,— the spin configuration at low temperatures 
can be determined by considering this shift in the light of 
Eqs. (l)-(7). We describe this procedure in more detail, 
which leads us to conclude the field-induced non-coUinear 
state, one of the findings in the present study. Since the 
applied field of 7.4 T (||a) is far greater than the spin-flop 
field, the local moments lie nearly in the be plane with 
a slight canting toward the field direction (||d); more- 
over, the AF interaction between Al and A2 makes the 
b and c components of local moments antiparallel; thus 



-<2 



and 



Ml, ^ 



^M 



A2- 



From these con- 



straints and Eq. (6), it turns out that the local moments 
cannot have the b component, i.e., M^, = M^, = 0. 
Thus [see Eq. (1)] is given by 



T^a A aa 



Ml, 
H 



j^ac Ml, 



^AV 



H 



Al- 



(8) 



While the sign of magnetization component parallel to 
the field. Ml,, is obviously positive, that of the per- 
pendicular component. Ml,, is not trivial. However, it 
can be judged explicitly from the sign of the shift at low 
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FIG. 4: (Color online) (a) Experimental "C-NMR shift of the 
outer site under 7.4 T parallel to the a axis, (b) Spin structure 
at low temperatures deduced from (a), (c) Hypothetical ^^C- 
NMR shift of the outer site and (d) spin structure at low 
temperatures deduced from (c). The horizontal lines in (a) 
and (c) represent the shift origin. In (b) and (d), the canting 
of local moments toward the a axis is exaggerated. 
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FIG. 5: (Color online) Spin susceptibility measured under 7 T 
and ^''C-NMR shift of the outer site. Solid lines are guide for 
the eyes. Inset: K-x plot. The hyperfine coupling constant 
was estimated from the linear fitting shown in the inset. 



temperatures, where |Mj^j^| 3> M'^^ > and therefore 
the second term is dominant in Eq. (8). Because A"^^ > 
as mentioned above (Sec. HI B), the large positive shift 
below Tn is attributed to the growth of positive M^^ . In 
this way, all components of the local moment at dimer Al 
at low temperatures are determined to be ^ 
> and = 0. By applying Eqs. (5)-(7) to Mai, 

one finally obtains the whole spin configuration at low 
temperatures, which is roughly sketched in Fig. (H^b). 

Figure Ufa) represents an anomalous feature over- 
looked in the previous studies: on cooling, the NMR 
spectra always show the positive shift, i.e., the negative 
shift as shown in Fig. HJc) has never been observed. 
This indicates that at low temperatures the spin config- 
uration having the negative M^^-^ as shown in Fig. IH^d) 
never appears. Note that such situation is not expected 
in the conventional AF transition under a magnetic field, 
where the Figs. HJ^b) andSfd) configurations are degen- 
erate and thus have the equal chance to emerge below 
Tn- Therefore, the present result implies that the for- 
mer configuration has a lower energy than the latter. As 
seen in Sec. V, the origin of this energy difference is 
actually understood by considering the DM interaction 
with the antisymmetric nature, i.e. D12 • {S1XS2) ^ 
D12 • (S'2x5'i). The point is that the thermal average 
of the two configurations yields positive M^^ (staggered 
component) as well as positive M'^^ regardless of tem- 
perature because of the energy difference, whatever the 
origin is. In this way, the inevitable positive shift in Fig. 



IHa) turns out to indicate that the field-induced local mo- 
ments under H\\a are not uniform (Mf > and M[ = 
0) but non-coUinear (Mf > and Af^^ = "M^^ > 0) 
even above zero-field Tn. 

The non-collincar spin state can be identified more evi- 
dently when the uniform Xspin is compared with the NMR 
shift K'^^ . Figure [5] shows the temperature dependence 
of NMR shift (if||a) and in-plane Xspin- We assume that 
the in-plane anisotropy of Xspin is negligibly small, as 
is the case in most of the organic conductors. There- 
fore, the uniform Xspin can be interpreted as probing 
M^j^. Note that when the field- induced local moments 
under H\\a are uniform (M° > and = 0) as usual, 
Xspin and K"^^ should show the same temperature de- 
pendence, because is given by A^^M^^/T? -I- C^i 
in this case [see Eq. (8)]. As seen in Fig. [5l however, the 
two quantities exhibit quite different temperature depen- 
dence above the zero- field Tn. Especially, the difference 
is pronounced in 30-60 K: the shift increases on cool- 
ing, indicating the growth of local magnetization, while 
the uniform Xspin (i-e., -^^i) decreases. This remarkable 
disagreement demonstrates that the finite staggered mo- 
ment M^j^ (= > 0) contributes to the observed 
shift in 30-60 K. 

Using the reported data of A'X-^ = 1.26^^° and C^i = 
96 (ppm) ■ — 1 — we can estimate the magnitude of stag- 
gered moment semi-quantitatively as follows. Here, Kj^^ 
[ppm] is given by (M^^/i? + 1.26M%^/H)+ 96. At 
high temperatures, it is plausible that the field- induced 
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FIG. 6; (Color online) Estimated temperature dependence of 
(a) the uniform and staggered moment under 7.4 T and (b) 
the non-coUinear angle 9 measured from the magnetic-field 
axis. Eye guides are shown together. 



local moment is nearly parallel to the magnetic-field (i.e., 
M^-^ » -^Ai)- Therefore Kj^^ around room temperature 
is given by A'^^M%^/H+ 96, and A"^^ can be estimated 
roughly from the K-x plot of the high-temperature data, 
as shown in the inset of Fig. [5l Note again that the 
difference between the observed shift and ^^"M^^/if-l- 
96 (i.e., the straight line in the inset) is attributed to 
the contribution of staggered component, A'^^M'^^ (= 
1.26A^"^Mj^j); thus, using the obtained one can 

evaluate the temperature dependence of the staggered 
component M"^^ as well as that of the uniform compo- 
nent M^^. The results are shown in Fig. [DJa), where 
the staggered component Mj^^ (= —M'^2) starts to grow 
markedly below ~ 60 K. The temperature dependence 
of non-coUinear angle 9 defined by arctan(Mj^j^/M^-^) is 
shown in Fig. El^b), which demonstrates that with low- 
ering temperature the field-induced non-coUinear local 
moments gradually tilt in accordance with the growth of 
staggered component. 

Note again that the positive staggered moment is in- 
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FIG. 7: (Color online) (a) Schematic illustration of the 
present square-lattice model in the presence of uniform and 
staggered fields. The local moment at site interacts with 
four nearest-neighbor sites, 1-4. -I- and — denote the two mag- 
netic sublattices (see the text), (b) Definition of the uniform 
moment Af(0), staggered moment M{Q) and non-coUinear 
angle 9. 



duced by a magnetic field even above "Tn". Therefore 
the AF symmetry breaking occurs in the whole tem- 
perature range under a field, although its degree is ex- 
pected to vanish asymptotically toward high tempera- 
tures. Thus, the rapid increase in the staggered moment, 
which has so far been addressed nominally as the AF 
transition, is not a phase transition but a crossover with- 
out thermodynamic singularity. As seen in Sec. V, the 
crossover behavior under a magnetic field is also demon- 
strated in the molecular-field calculations. 



V. NUMERICAL RESULTS 
— MOLECULAR-FIELD APPROXIMATION 

A. Model Hamiltonian 

To reveal the origin of the unconventional magnetism 
highlighted above, we investigate the S = 1/2 AF Heisen- 
berg model including the DM and Zeeman interactions 
on a two-dimensional square lattice [Fig. El^a)]. Here, 



1-1 = J ^ Si - Sj+ gfiB ^Si- 

<ij> i 

+ J2 ■ {S. X S,l (9) 

<ij> 

where Si and Hi are the spin operator and external mag- 
netic field, respectively, at site i, and < ... > denotes the 
sum over the nearest-neighbors. For simplicity, the DM 
vector D is set to be parallel to the z axis, and we ignore 
the AF exchange between the next-nearest-neighbors. A 
uniform field H{q = 0) is applied along the y axis, where 
q is a wave-number vector. In the case that only the AF 



exchange term is present (i.e.. Hi 



D, 



0), the Neel 



order with Q = (tt, tt) is emergent at low temperatures; 
thus the square lattice can be divided into two magnetic 
sublattices, -I- and — , as shown in Fig. [Tfa). We choose 
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FIG. 8: (Color online) Molecular-field calculation of temperature dependence of (a) the uniform magnetization M(0) and the 
staggered magnetization M{Q) under 15 T, (b) the staggered susceptibility xiQ) under OT and 15 T, and (c) the non-coUinear 
angle 6 under 15T. In the present calculations, the DM interaction is not included. 



the X axis as the easy axis, but the magnetic anisotropy 
is not taken into account explicitly. 

In this section, we solve the present model using the 
molecular-field approximation on the two magnetic sub- 
lattices. The molecular-field Hamiltonian of Eq. (9) is 
given by 



cfF 



5MB 



■off 



(10) 



Here i+ {i—} denotes a lattice point on the -I- (— ) mag- 
netic sublattice, and 11°^ {H°_^) is the effective field per- 
ceived by the + (— ) magnetic sublattice. The effective 
field is described by the sum of the external and molec- 
ular fields: 



7 T 7 

Hf = -^^(M_) - ^(M_) X D+ + H,, 



7 J 7 

H^^ = -^^{M+) - ^£>+_ X (M+) + H,. 



(11) 

i • 

(12) 



Here is along the z axis, Z = 4 is the number of the 

nearest- neighbors, and ((TW-)) denotes the statis- 

tically averaged local moments on the + (— ) sublattices. 
In the molecular-field solutions, and (M^) are 

determined self-consistently so as to satisfy {M-^-)\\H'^ 
and {M^)\\H''_^. Note that because of Zeeman and DM 
terms such local moments are symmetric with respect to 
the y axis [Ijiii'(O)] and lie in the x-y plane. For clarity, we 
adopt the basis composed of the uniform moment along 
the y-axis, M(0) = |(M+) + (M_)|/2, and the staggered 
moment along the x-axis, M{Q) = |(M+) - (M_)|/2, 
as illustrated in Fig. [T^b). We define the non-collinear 
angle 9 as arctan~-'^[M(Q)/M(0)]. When we calculate 
the staggered susceptibility x{Q) under H{0), a stag- 
gered field H{Q) is applied additionally along the x axis 
[see Fig. El^a)]; then x{Q) is obtained from the difference 
of M{Q) between under H{Q) and under H{0)+H{Q), 



as represented in the following equation: 

_ M{Q,H{0) + H{Qj)-M{Q,H{0)) 

Although x{Q) should be calculated under an infinites- 
imal H{Q), the present calculation of x{Q) was per- 
formed under \H{Q)\ = 20 mT because of the finite 
numerical accuracy. Below, within the molecular field 
approximation, we show the temperature dependence 
of M{0), M{Q) and x{Q) under various magnitudes 

of HiQ) and Throughout the calculation, J 

20 K and g — 2 were used. The molecular-field Tn 
in the square lattice under a zero field is given by 
y/.P + I-D+-P. Hence in the case of = 0, Tn 

equals J, namely, 20 K. 



B. Heisenberg model without DM interaction 
under magnetic field 

It is instructive to show first the consequences of the 
Heisenberg model without the DM interaction under a 
magnetic field. As is well known, the AF transition, 
which is characterized by the emergence of M{Q) and 
the divergence of x{Q)i is still well-defined even under 
a uniform field of moderate magnitude. This feature is 
reproduced in our numerical calculation: under 15 T, fi- 
nite M{Q) appears below Tn (- 19.57 K) [Fig. E^a.)], 
and xiQ) is divergent at this temperature [Fig. Eljb)], 
although Tn is suppressed slightly from 20 K.^^ Below 
Tn, M{Q) can either be positive or negative (i.e., par- 
allel or antiparallel to the x axis), because the positive- 
and the negative-M(Q) phases are degenerate. The pres- 
ence of this spontaneous symmetry breaking also assures 
the well-defined AF transition even under 15 T. Reflect- 
ing the well-defined AF transition, the non-collinear state 
(i.e., finite 9) is emergent only below Tn, as shown in Fig. 
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FIG. 9: (Color online) Schematic configuration of the local 
moments, molecular-fields produced by them, and external 
uniform field for the case of (a) coUinear local moments and 
(b) non-colhnear local moments. The effective fields H"^ 
(if?) perceived by the + (— ) sublattices are also shown for 
each case. 



C. Heisenberg model with DM interaction under 
magnetic field 




FIG. 10: (Color online) Molecular-field calculations of tem- 
perature dependence of (a) the staggered magnetization 
M{Q) and (b) the non-colhnear angle 9. Insets: Plots in 
(a) linear scale and (b) logarithmic scale. 
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Next, we consider the Heisenberg model including the 
DM interaction under a uniform field H{0) within the 
molecular-field approximation. The model was found to 
exhibit the following characteristic features: (i) a uniform 
field induces the non-colhnear local moments even in the 
paramagnetic phase (i.e., above the zero-field In); (u) the 
divergence of x{Q) at Tn is suppressed under a uniform 
field; and (iii) the peak temperature of xiQ) increases 
with the field in a low-field region. Below we explain 
these properties in sequence. 

Under a zero field, the self-consistent moments are 
zero at high temperatures (i.e., paramagnetic), while the 
spontaneous local moments appear below the mean-field 

Tn (= + In the paramagnetic state, the 

local moments are induced under a finite field. Note that 
such field-induced local moments cannot be coUinear, be- 
cause the coUinear state does not satisfy {M-^-)\\H'^ and 
(M_)||i?iff, as illustrated in Fig. H^a). Therefore, in 
a self-consistent result, the field-induced local moments 
should form a non-coUinear configuration characterized 
by positive M(0) and positive M(Q)[Fig. M.^)] in the 
whole temperatures. We actually calculated the tem- 



perature dependence of the field-induced local moments. 
Figures fTUf a) andfTUTb) show the behavior of M{Q) and 
6*, respectively, under |i?(0)| = 3 T with = 0.2, 

1, and 3 K. Although M{Q) and 9 show rapid increase 
around 20 K on cooling, they are finite in the whole 
temperature range as expected, demonstrating that the 
AF transition accompanied by thermodynamic singular- 
ity does not occur in the presence of DM and Zeeman 

interactions. We also note that as |il>-| | increases, the 

non-coUinear feature above ~ 20 K becomes more pro- 
nounced over a wide temperature range [Fig. fWb)].^^ 

Figures fTlTa) and fTlTb) show the temperature depen- 
dence of x{Q) for I = 0.2 and 3 K, respectively, 

under various H{0). As H{0) increases, the xiQY 
divergence emergent at T is gradually suppressed into a 
rounded peak, consistent with the absence of AF transi- 
tion. The peak suppression is more prominent for larger 

I, as seen in Fig. [TIT b). These results arc in sharp 

contrast to the magnetic-field effect on x{Q) without the 
DM interaction [Fig. W)^)]- 

From the experimental point of view, the peak struc- 
ture of x{Q) [or x(Q)-related quantities] has often been 
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FIG. 11: (Color online) Molecular-field calculations of tem- 
perature dependence of x(Q) under various magnetic-fields at 
(a) \D+- \ = 0.2 and (b) \D+- \ = 3 K. Inset of (b): Enlarged 
view of low-x(<5) part. 



interpreted as the indication of AF transition. How- 
ever, we emphasize again that in the presence of DM 
interaction the peak under a magnetic field should be 
regarded as the manifestation of the crossover. There- 
fore we name the peak temperature the crossover tem- 
perature, T*. The magnetic-field dependence of Tn (at 
\D+_\ = 0) and T* (at ^ 0) is quahtatively dif- 

ferent, as shown in Figs. HDJa) and [T^ b): T* increases 
with H{0) in a low-field region regardless of the strength 
of , while Tn decreases monotonously. The field- 
dependence of T* is determined by the superposition of 
the following two effects. One is the competition between 
AF exchange and Zeeman interactions, which leads to 
the suppression of Tn, as is well known. The other is 
the cooperation between the DM and Zeeman interac- 
tions, which is naturally expected from the fact that, 
compared with a coUinear-AF state, a canted- AF state is 
favored by both interactions. This cooperation effect ap- 
pears in AT = T* — Tn, which is plotted against \H{0)\ 
in Fig. [T^ c). In contrast to the case of Tn, AT increases 
monotonously with \H{0)\, and this increase is enhanced 
for larger | value. The field dependence of T* re- 



20.2 




H(T) 



FIG. 12: Molecular-field calculations of magnetic-field depen- 
dence of r* at (a) = 0.2 and (b) j£>+-| = 3 K. (c) 
Molecular-field calculations of magnetic-field dependence of 
AT (= T* - Tn) at \D+-\ = 0.2 and 3 K. For comparison, 
the magnetic-field dependence of Tn at = K is also 
shown in (a) and (b). 

suits from the sum of the two effects and thus may show 
non-monotonic field dependence such as shown in Fig. 

VI. DISCUSSION 

A. Comparison between the experimental and 
numerical results 

Now we compare the experimental results with 
molecular-field calculations. The puzzling magnetism 
found in k-CI is summarized as follows: (i) even above 
the zero-field Tn , the field- induced local moments have a 
staggered component perpendicular to the applied field, 
and thus the AF transition is not well defined under a 
magnetic field; (ii) with increasing magnetic field, the 
peak structure in the temperature dependence of 1 /Ti is 
suppressed;^^ and (iii) the peak temperature in 1/Ti in- 
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creases with the apphed field.— As seen in Sec. V, the ex- 
perimental result (i) was reproduced qualitatively by the 
present calculation. In addition, the calculations showed 
that the peak structure of xiQ) is suppressed with in- 
creasing I J? (0)1 (Fig. [TT|) and that the peak temperature 
increases in a low-field region [Figs. fT2ya) and fT2yb)]. 
Thus, as far as l/TiT around Tn (or T*) reflects well 
x(Q),^"' the molecular-field results also explain features 
(ii) and (iii) qualitatively. These agreements between the 
experimental and molecular-field results indicate that the 
puzzling magnetism in k-CI is essentially due to the in- 
terplay between Zeeman and DM interactions. It should 
be noted that the NMR line shift in the spin system in- 
cluding the DM interaction is not proportional to the 
macroscopic magnetization even in the nominally param- 
agnetic state except the high-temperature regime; as was 
demonstrated by our experimental and numerical results, 
a magnetic field induces the non-coUincar spin configu- 
ration and thus the NMR line shift reflects not only the 
uniform component but also the staggered component of 
local moments. 



order and thus x('Zaf) generally grow toward Tn. This 
growth may enhance the effective magnitude of the DM 
interaction, resulting in a more prominent staggered mo- 
ment and in more field-dependent T* than expected in 
the molecular-field calculation. 

The molecular-field results do not reproduce the mono- 
tonic decrease in Xspin down to ~ 30 K in k-C1 [the inset 
of Fig. [5] even in a qualitative level. This discrepancy 
is possibly due to the low dimensionality of k-C1, and 
the further investigation of Eq. (9) is needed beyond the 
molecular-field approximation. On the other hand, it is 
likely that the magnetism of k-C\ can not be described 
adequately by the model Hamiltonian of Eq. (9), be- 
cause K-Cl is located near the Mott transitionj^ii^i^iiSiiii 
thus the charge fluctuations are probably not negligible 
above the critical endpoint at ^ 40 K. Therefore the lo- 
calized spin model may not be adequate to describe the 
magnetism of k-C1 at high temperatures. The Hubbard 
model, which incorporates the charge degrees of freedom, 
or the Heisenberg model including higher-order exchange 
interactions are considered more realistic. 



B. Quantitative disagreement between the 
experimental and molecular-field results 

The agreement between the unconventional magnetism 
of K-Cl and the molecular-field results is qualitative, but 
not quantitative; the features revealed by the calcula- 
tions are highly enhanced in the experimental observa- 
tions. In K-Cl, D/J is considered to be 10"'^-10~^ (Ref. 

19); therefore the molecular-field results for \D^ | — 0.2 

K (hence D/J= 10~^) should be a relevant reference. 
However, the unconventional features in the numerical 

results for |£>_| | = 0.2 K are not as pronounced as those 

observed in k-C\. First, the values of 6 above Tn in the 
numerical results [Fig. [TUTb)] are too small to be ob- 
served experimentally, while the experimental 9 values 
remain observable even well above Tn [Fig. E^b)]. Sec- 
ond, the increase in T* with a uniform field is also tiny 
(0.7 % at most) in the calculations [Fig. [T^ a)]. while 
an increase by more than 20 % was observed in k-C1.^* 
When these quantitative deviations are considered from 
the molecular- field point of view, D/J of k-C1 seems as 
if it were much larger than 10~^. These quantitative dis- 
agreements imply that some sort of effects beyond the 
molecular-field treatment enhances the unconventional 
features induced by the interplay between the DM and 
Zeeman interactions. Although we have no numerical ev- 
idence, we speculate that the most possible candidate is 
the AF short-range order enhanced by low dimension- 
ality, considering that k-C1 is actually a quasi-2D sys- 
tem. In fact, the values of Tn/J differ significantly for 
the molecular- field and experimental results: Tn/J = 1 
in the present molecular-field results, while Tn/J is less 
than 0.1 in k-C1 (the zero-field Tn is '-^ 23 K and J is 
considered to be of the order of room temperature) 
In low-dimensional AF spin systems, the AF short-range 



C. Comparison with the previous studies 

Our explanation of the field-dependent 1/Ti in k-C1 
is different from the one proposed by Hamad et al.^^ In 
general, l/TiT is given by 

TiT ^ u;o ^ ^ 

where loq is an observing NMR frequency proportional 
to |i?(0)|. In most cases, Im x(q',ix>)/u; is w-independent 
at low frequencies of NMR and hence 1 /Ti does not de- 
pend on a magnetic field. In order to explain the strongly 
field-dependent 1/Ti of k-CI, Hamad et al. argued that 
Im x(q, uj)/u! can have w-dependence due to slow spin dy- 
namics under spin frustration. This is a frequency effect, 
but not a field effect (i.e., the magnetic-field dependence 
of x(Q,t<^=0) is not considered). They performed numer- 
ical calculations on the Heisenberg model (without the 
DM interaction) on an anisotropic triangular lattice by 
using the modified spin- wave theory and reproduced the 
suppression of 1/Ti peak under a uniform field. However 
their numerical results also exhibit three characteristics 
inconsistent with the experiments as follows: a double- 
peak structure appears in the temperature dependence 
of 1/Ti; the peak temperature of 1/Ti does not depend 
on |Jf(0)|; and 1/Ti depends on a magnetic field over a 
wider temperature range than observed. 

In contrast, we ignore the frequency effect and consider 
the field dependence of x(q, uj—0) in the presence of DM 
interaction. This is the field effect. Our treatment ex- 
hibits the strong field dependence of x{Q) only around 
the zero-field Tn [Figs. [TlTa') andfTlTb')]. As described in 
Sec. VI A, as far as x{Qaf) 1/TiT are well correlated 
near Tn (or T*) in k-CI^^ the molecular-field results ex- 
plain the experimental results without the discrepancies 
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encountered in the model postulated by Hamad et alJ^ 
Thus, the unconventional field dependence of 1/Ti ob- 
served only around the zero-field Tn is successfully un- 
derstood in terms of the field effect on x(9af)- 

The field-induced staggered moment, which is a key 
consequence of the present study, has often been observed 
and discussed in quantum disordered or gapped spin sys- 
tems. Oshikawa and Affleck'^^ argued that in S = 1/2 AF 
spin chains, the interplay between the Zeeman and DM 
interactions induces the staggered field perpendicular to 
the applied magnetic field, resulting in the field-induced 
staggered moment. This actually explains the observa- 
tions in Cu benzoate^ and copper pyrimidine dinitratei^ 
The spin-gapped system with DM interaction, such as 
NENP (S* = 1 AF spin chain)^>29 and SrCua (603)2 
{S = 1/2 2D frustrated spin system) also show a 
field-induced staggered moment below a critical field. 
These properties are also well understood in terms of 
the staggered-field mechanism. The present study is ad- 
dressed as the extension of this issue to the quasi-2D spin 
systems undergoing a finite-temperature AF transition. 
Recently, the NMR line shift of La2Cu04, in which the 
DM interaction is inherent, has been discussed in terms 
of the field-induced staggered moment/^ 

VII. CONCLUSION 

We performed the magnetization and NMR measure- 
ments on K-(BEDT-TTF)2Cu[N(CN)2]Cl for deeper un- 



derstanding of its magnetism. We found that the ap- 
plication of a uniform field induces a staggered moment 
perpendicular to the applied field even above Tn. The 
molecular-field calculation revealed that the field-induced 
staggered moment comes from the interplay between the 
Zeeman and DM interactions and that this interplay ex- 
plains the several otherwise puzzling magnetic features 
observed in experiments. By comparing the experiment 
results with our numerical results, we concluded that 
the DM interaction affects significantly the field-induced 
magnetism in a paramagnetic phase especially in low- 
dimensional systems. 
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